DiscreteLyapunovSolve
DiscreteLyapunovSolve[a,c] finds the numeric solution x of the discrete matrix equation a.x.a-xc.
DiscreteLyapunovSolve[a,b,c] solves a.x.b-xc.
DiscreteLyapunovSolve[{a,d},c] solves a.x.a-d.x.dc.
DiscreteLyapunovSolve[{a,d},{b,e},c] solves a.x.b-d.x.ec.
Please visit the official Wolfram Language Reference for more details and examples on core symbols.